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Collective behavior of "electronic fireflies" 
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A simple system composed of electronic oscillators capable of emitting and detecting light-pulses 
is studied. The oscillators are biologically inspired, their behavior is designed for keeping a desired 
light intensity, W , in the system. From another perspective, the system behaves like modified 
integrate and fire type neurons that are pulse-coupled with inhibitory type interactions: the firing 
of one oscillator delays the firing of all the others. Experimental and computational studies reveal 
that although no driving force favoring synchronization is considered, for a given interval of W 
phase-locking appears. This weak synchronization is sometimes accompanied by complex dynamical 
patterns in the flashing sequence of the oscillators. 
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I. INTRODUCTION 

Synchronization of quasi-identical coupled oscillators is 
one of the oldest and most fascinating problems in physics 
[ESQ- Its history goes back to C. Huygens who first no- 
ticed the synchronization of pendulum clocks hanging on 
the same wall. Besides mechanical or electric oscillator 
systems, nature is also full with several amazing examples 
in this sense d,[5, 6]. Synchronization in all these systems 
appears as a result of some specific coupling between the 
units. This coupling can be local or global, and can be re- 
alized through a phase-difference minimizing force @, [1] 
or through the pulses emitted and detected by the os- 
cillators [3, [13] • In most of these synchronizing systems 
there is a clear driving force favoring synchronization, 
and in such way the appearance of this collective behav- 
ior is somehow trivial. In the present work however, a 
nontrivial synchronization will be presented. This weak 
synchronization (phase-locking) appears as a co-product 
of a simple collective optimization rule. 

One well-known phenomena which inspired us in this 
work is the collective behavior and synchronization of 
fireflies [llj. Although our aim here is not to model fire- 
flies, the oscillators ("electronic fireflies") considered in 
our system are somehow similar to them: they are capa- 
ble of emitting light-pulses and detecting the light-pulse 
of the others. In this sense our system is similar to an en- 
semble of fireflies although the coupling between the units 
is different. From another perspective, the oscillators be- 
have like pulse-coupled " integrate and fire" type neurons 
0) OH Q3- Contrary to the classical integrate and fire 
oscillators, in the considered system an inhibitory type 
global interaction is considered. This means that the fir- 
ing of one oscillator delays (and not advances) the phase 
of all the others. This system does not necessarily fa- 
vor synchronization, it is rather designed to keep a de- 
sired W light intensity in the system. This light intensity 
is controlled by a firing threshold parameter G imposed 
globally on the oscillators. Surprisingly, as a co-product 
of this simple rule, for certain region of the firing thresh- 
old parameter phase-locking and complex patterns in the 




FIG. 1: Experimental setup. The photo on the left shows the 
"electronic fireflies" (oscillators) placed on the circuit board. 
The photo on the right shows one oscillator. 



flashing sequence of the oscillators will appear. We be- 
lieve that such dynamical laws could be realistic for many 
biological systems. 

The studied system will be described in more details in 
the following section. The used electronic device will be 
briefly presented and the obtained non-trivial collective 
behavior will be studied. In order to get more confidence 
in the observed non-trivial results computer simulations 
were also performed. 



II. THE EXPERIMENTAL SETUP 

The constructed units are integrate and fire type oscil- 
lators Q with a modified interaction rule. Their coupling 
and communication is through light, the units are capa- 
ble of emitting and detecting light-pulses. The oscillators 
are practically realized by a relatively simple circuit, the 
main active elements being a photoresistor and a Light 
Emitting Diode (LED). Each oscillator, i, has a charac- 
teristic voltage Ui, which depends on the resistance, Ri, 
of its photoresistor. The global light intensity influences 
the value of Ri in the following sense: when the light in- 
tensity increases Ri decreases, leading to a decrease in Ui. 
In the system there is a global controllable parameter G, 
identical for all oscillators. By changing the parameter 
G, one can control the average light intensity output, W, 
of the whole system. If the voltage of the oscillator grows 
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FIG. 2: Circuit diagram of one oscillator. 



above this threshold (Ui > G) the oscillator will fire, this 
meaning its LED will flash. This flash occurs only if a 
minimal time period T m i ni has expired since the last fir- 
ing. The oscillator has also a maximal period, meaning if 
no flash occurred in time T maXi , then the oscillator will 
surely fire, fn laymen terms firing is favored by darkness 
and the value of the controllable G parameter charac- 
terizes the "darkness level" at which firing should occur. 
Through this simple rule the G parameter controls the 
average light intensity output of the system. The tech- 
nical realization of the above dynamics is illustrated in 
Fig. O After the system has fired the 22 /iF capaci- 
tor is completely discharged by the negative pulse from 
the inverted output of the monostable. As soon as the 
light flash ended, the same capacitor will start charging 
from the current flow through the 270 KS1 resistor. The 
IC1B comparator will trigger another flash as soon as 
the potential on the mentioned capacitor will overcome 
the value fixed by the group of three resistors on its pos- 
itive input (the firing threshold). Two of the resistors 
are connected to constant potentials (ground and +5 V), 
the third resistor is connected to the output of the sec- 
ond comparator IC1A which will have a value depending 
on the ratio between the reference voltage and a certain 
amount of light measured by the photo resistor. The flash 
time is determined by the second capacitor together with 
the 12 Kf2 resistor connected to the monostable. The 
photoresistor has a relatively low reaction time around 
40 ms, while the minimal and maximal period of firing 
are around 800 ms and 2700 ms. The time of one flash 
is around 200 ms. 

The oscillators are placed on a circuit board in the form 
of a square lattice (see Figured]). The maximal number 
of oscillators which can be included are 24. A computer 
interface and program controls the G threshold parame- 
ter and allows us to get information automatically about 
the states of all oscillators. The state of an oscillator 
is recorded as if the oscillator does not emit light and 
1 when the oscillator fires (emits light). Whenever the 
state of the oscillator system changes, the program writes 
in a file the corresponding time with a precision of mil- 
liseconds and the new states of the units. 

To obtain an enhanced global interaction the whole 
system is placed inside a closed box. The box has mat 



glass mirror walls to uniformly disperse the light-pulses in 
the box. A graphical interface allows to visually control 
the state of the units. 

In order to fully understand the behavior of the sys- 
tem one has to accept that the coupling between pairs of 
oscillators are not exactly of the same strength. Also, 
the characteristic electronic parameters differ slightly 
(2 — 10%) among the units. 



III. COLLECTIVE BEHAVIOR 

At constant light intensity one unit behaves as a simple 
stochastic oscillator. Whenever the G threshold is under 
a given G c value the oscillator will fire with its mini- 
mal period and above G c with its maximal period. G c 
depends of course on the imposed light intensity. Con- 
sidering more oscillators (i = 1, ...,n) and by letting 
them interact, interesting collective behavior appears for 
a certain range of the G threshold parameter. 

Due to the inhibitory nature of the considered inter- 
action, during the firing of oscillator i the characteristic 
voltages of the others (Uj, j ^ i) will decrease. If the G 
parameter is so small, that under this condition the other 
oscillators can still fire (Uj > G), than all oscillators will 
fire in an uncorrelated manner. Each of them will be fir- 
ing at its own T m i n . period and the interaction is thus 
not efficient. In such case no collective behavior can be 
observed. 

Increasing the value of G will make the pulse- like inter- 
action efficient. The oscillators will avoid firing simulta- 
neously and a simple phase- locking phenomenon appears. 
The pulse of one unit (let us assume i) delays the fir- 
ing of the others by decreasing their voltages below the 
threshold: Uj < G,j ^ i. Due to the tiny differences 
in the coupling between the pairs (caused for example 
by different distances) and in the parameters of the elec- 
tronic elements, the Uj voltages are different. The imme- 
diate consequence of this is that the next firing will occur 
most probably in the oscillator with the highest voltage 
(counting of course only those oscillators, which are al- 
ready capable of firing). This oscillator is the one which 
was influenced the less by the light-pulse of the previ- 
ous firing. If the total combined time of firing for the 
n oscillators is smaller than the period T max the result 
is that after very short time phase-locking appears and 
a firing chain (with period T 6 [T m i n , T max ]) will form, 
each oscillator firing in a well-defined sequence. For a 
given system and a fixed G threshold this stable pattern 
is always the same. If the total time of firing of the n 
oscillators exceeds T maxi the firing pattern will be much 
longer and more complex. 

Increasing further and over a limit the G threshold pa- 
rameter the previously discussed weak synchronization 
(phase-locking) disappears. In this case the voltages of 
all oscillators are much smaller than the threshold value 
Ui < G, so the firing of a unit can not influence the oth- 
ers. All oscillators will fire with their own T maXi period 
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and no interesting collective behavior is observed. Again, 
the interaction is not efficient. 

The collective behavior of the system can be easily 
analyzed by plotting a kind of phase-histogram for the 
oscillator ensemble. Choosing a reference oscillator, the 
relative phases of all the others are defined by measuring 
the time difference between their pulse and the last pulse 
emitted by the reference oscillator. Studying these time- 
delays during a longer time period a histogram is con- 
structed for their distribution. This histogram shows how 
frequently a given time-delay occurred and gives thus a 
hint whether a constant firing pattern is formed or not. 

Experimental and computer simulated results for the 
phase-histogram confirm the above presented scenario of 
the collective behavior. As an example, on Fig. [3l re- 
sults obtained on a relatively small system with n = 5 
oscillators are shown. In the first column of Fig. [3] (fig- 
ures a, b, c and d), experimental results for four differ- 
ent values of the G threshold are plotted. For a small 
threshold parameter (G = 500 mV), no self-organization 
appears (Fig. (3^)- Due to the fact that the character- 
istic time-periods of the oscillators are slightly different, 
almost all values will occur with the same probability 
in the phase-histogram. Beginning with G — 1300 mV 
a kind of order begins to emerge, and a trend towards 
the self-organization of the oscillator pulses is observed 
(e.g. Fig. [3b for G = 2000 mV). In the neighborhood 
of G = 3000 mV threshold value (Fig. [3]:) clear phase- 
locking appears. One can observe that a stable firing 
pattern has formed, each oscillator has an almost exact 
phase relative to the reference oscillator. For an even 
higher value (e.g. G — 4200 mV), disorder sets in again, 
phase-locking disappears and all oscillators fire indepen- 
dently with their own maximal period (Fig. [3Ji). 

In the second column of Fig. [3J we present the corre- 
sponding simulation results. In simulations the parame- 
ters of the oscillators are defined as following: the average 
minimal time period is T m i ni = 900 ms, the average max- 
imal period T maXi = 2700 ms, and the average flashing 
time T fi as h — 200 ms. For an easier comparison, the val- 
ues are chosen to be similar with the real experimental 
data. We considered a uniform distribution of the oscilla- 
tors parameter around these average values using a ±50 
ms interval for T m i n and T max and a ±20 ms interval 
for T fi as h- One could argue of course that a Gaussian 
distribution would be much more appropriate, but given 
the fact that we simulate here relatively small number of 
oscillators the exact statistics is irrelevant. Considering 
some deviations from the average is however important 
in order to reproduce the collective behavior of the sys- 
tem. An uncorrelated noise in time is also considered. 
This will randomly shift the T min , T max and Tfi as h pe- 
riods of each oscillator at each cycle. Again, a uniform 
distribution on a ±20 ms interval was considered. The 
characteristic voltages of the oscillators are set to be in 
the interval 4100 ± 100 mV in dark, 2100 ± 100 mV when 
one single LED is flashing and 1050 ± 100 mV when two 
LEDs are flashing simultaneously. Whenever k LEDs are 
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FIG. 3: Relative phase histogram for n = 5 oscillators. Exper- 
imental results are in the first column, and the corresponding 
simulation results are in the second column. Figures a) and 
e) are for G = 500 mV; figure b) and f) are for G = 2000 mV; 
figures c) and g) are for G = 3000 mV, and figures d) and h) 
are for G = 4200 mV. 



simultaneously flashing the characteristic voltages of the 
others are considered to be 2100/A:±100 mV, however for 
n = 5 oscillators only very rarely happens to have more 
than two oscillators simultaneously firing. This values 
were chosen to approximately match the experimental 
ones, and we do not try here to give a theoretical model 
for the nonlinear behavior of the photoresistor. Fluctua- 
tions in time and among the parameters of the oscillators 
are again included. Differences in the strength of the cou- 
pling between pairs of oscillators are however neglected. 
Using these parameters, it is assumed that each oscilla- 
tor can flash whenever its voltage exceeds the threshold 
G. The flashing cannot occur earlier than T m i ni or later 
than T maXi relatively to its last firing. On Fig. [3^.,f.,g. 
and h., the simulated phase- histograms of the oscillators 
are plotted and compared with the corresponding exper- 
imental data. The observed experimental results, includ- 
ing the non-trivial synchronization (phase- locking) , were 
successfully reproduced. 

It is also possible to define a kind of order-parameter 
that characterizes the observed synchronization level. 
Our method for calculating this is the following: 

1) A reference oscillator k is chosen and the phases of 
all oscillators are calculated relative to this oscillator. 

2) Let hi{f) denote the value of the normalized phase- 
histogram for oscillator i (i = 1, . . . , n, i ^ k) correspond- 
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FIG. 4: Order parameters calculated from experimental (cir- 
cles) and simulation (dashed line) results plotted as a function 
of the G threshold. Systems with n — 3, 5, 7, 9 oscillators are 
considered. 

ing to phase difference (time difference) value /. Since 
we have a normalized histogram, hi(f) G [0, 1] gives the 
occurrence probability of phase difference value / during 
the measurement (J^f hi(f) — !)■ 

3) A window of width a is defined (we have chosen a = 30 
ms). Shifting the window with A/ = 1 ms step, for each 

discretized value of / the sum Hi(f) — Ej=/-o/2 
is calculated for each oscillator i. 

4) Let r/j denote the difference between the maximum 
and minimum value of H{(f) averaged over all oscilla- 
tors: r k = YJi=i,(i^k) max(Hi) - min(Hi). 

5) Items 1-4 are repeated considering each oscillator in 
the system as reference oscillator. 

Finally, an averaging is performed over all the obtained 
rk values (fc = 1, . . . , n). The final order parameter is 
calculated thus as r = (rfe)fe. Averaging as a function 
of the reference oscillator is beneficial in order to get a 
smoother curve when only partial phase locking is de- 
tected (Figure [3)3.). In such cases the phase-diagrams 
are very sensible on the choice of the reference oscillator. 
On Fig. [2] the r order parameter is plotted as a function 



of the G threshold value. Systems with n = 3, 5, 7 and 
9 oscillators are considered. Experimental (circles) and 
simulation results (dashed line) are again in good agree- 
ment. The figure also illustrates that for an intermediate 
G interval value phase-locking appears. This weak syn- 
chronization is better (r is bigger) when there are less 
units in the system. One obvious reason for this is that 
by increasing n the total time of firing of the oscillators 
will increase and slowly exceed the value T max . As a re- 
sult of this the firing pattern will change from a simple 
"firing chain" to a much longer and more complicated 
pattern, decreasing the value of the order parameter. 

From Fig. 2] it is also observable that the experimental 
results show more intensive fluctuations. The reason for 
this is probably the complex noise present in the system. 
IV. CONCLUSION 

A system of electronic oscillators communicating 
through light-pulses was studied. The units were de- 
signed to optimize the average light intensity of the emit- 
ted light-pulses, and no direct driving force favoring syn- 
chronization was considered. Although our experiments 
focused on relatively small systems (up to 24 oscillators) 
interesting and rich collective behavior was observed. As 
a nontrivial result it was found that the inhibitory cou- 
pling induced a partial phase-locking for a certain interval 
of the controllable threshold parameter. This weak syn- 
chronization was realized by complex flashing patterns 
of the units. We believe that this study inspires fur- 
ther interesting research projects in which separately pro- 
grammable oscillators will be studied with various inter- 
action rules. Many other interesting collective behaviors 
can be obtained by controlling individually the parame- 
ters of the units, their interaction rule and the type of 
coupling between them. Such systems could also yield 
a new approach to unconventional computing, being in 
many sense similar with the presently developed CNN 
computers [l3j . 
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